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ABSTRACT 

We report the results of XMM-Newton observations of HD 49798/RX J0648. 0-4418, the only known 
X-ray binary consisting of a hot sub-dwarf and a white dwarf. The white dwarf rotates very rapidly 
(P=13.2 s) and has a dynamically measured mass of 1.28±0.05 M@. Its X-ray emission consists of 
a strongly pulsed, soft component, well fit by a blackbody with YYbb ^40 eV, accounting for most 
of the luminosity, and a fainter hard power-law component (photon index *~1.6). A luminosity of 
~10 32 erg s _1 is produced by accretion onto the white dwarf of the helium-rich matter from the wind 
of the companion, which is one of the few hot sub-dwarfs showing evidence of mass-loss. A search 
for optical pulsations at the South African Astronomical Observatory 1.9-m telescope gave negative 
results. X-rays were detected also during the white dwarf eclipse. This emission, with luminosity 
2xl0 30 erg s _1 , can be attributed to HD 49798 and represents the first detection of a hot sub-dwarf 
star in the X-ray band. HD 49798/RX J0648. 0-4418 is a post-common envelope binary which most 
likely originated from a pair of stars with masses ~8-10 Mq. After the current He-burning phase, 
HD 49798 will expand and reach the Roche-lobe, causing a higher accretion rate onto the white dwarf 
which can reach the Chandrasekhar limit. Considering the fast spin of the white dwarf, this could 
lead to the formation of a millisecond pulsar. Alternatively, this system could be a Type la supernova 
progenitor with the appealing characteristic of a short time delay, being the descendent of relatively 
massive stars. 

Subject headings: binaries: close - subdwarfs: individual (HD 49798) - white dwarfs - X-rays: indi- 
vidual (RX J0648.0-4418) 



1. INTRODUCTION 

Interacting binaries containing accreting white dwarfs 
constitute a large fraction of the population of bright 
Galactic X-ray emitters. These systems, including both 
persistent and transient sources, display a rich variety 
of properties, determined mainly by the strength of the 
white dwarf's magnet ic field and the nature of the mass 
donor star (see, e.g., iKuulkers et al.l (|2006l )). Most ac- 
creting white dwarfs are found in cataclysmic variables, 
where accretion proceeds through Roche-lobe overflow. 
The mass donor in these systems is typically a main se- 
quence star or another white dwarf, but also sub-giant 
and giant companion stars have been observed. Sym- 
biotic systems, where a white dwarf accretes from the 
stellar wind of a massive supergiant, have also been de- 
tected as X-ray sources. 

The system discussed here is the only known X-ray bi- 
nary composed by a white dwarf (RX J0648. 0-4418) and 
a hot sub-dwarf star (HD 49798). Hot sub-dwarfs are 
evolved low mass stars that lost most of their hydrogen 
envelopes and are believed to be in the helium cor e burn- 
ing sta ge (for a recent review on hot sub-dwarfs see lHeberl 
(2009)). They are spectroscopically classified in different 
types: sdB (with effective temperature T~25,000-28,000 
K), sdOB (with T-33,000-40,000 K), and sdO (with 
T>40,000 K). One possible mechanism responsible for 
the loss of their hydrogen envelopes is non-conservati ve 
mass transfer in a binary (|Podsiadlowski fc Hanll200l . 



Being one of the brightest hot sub-dwarfs (apparent 
magnitude V=8.3), HD 49798 has been the object of 
several studies in the optical/UV wavebands. The first 
spectroscopic observations led to its classification as a 
sub-dwarf of 06 spectral type (fjaschek & Jaschckl fl963| ) 
and showed radial velocity variations, later found to be 
caused by orbital motion with a period of 1.55 days 
(|Thackeravlli970l iStickland fc Llovdl[l99l . No spectro- 
scopic evidence for a secondary could be found, leading 
to the early suggestion that the companion could be a 
white dwarf ()Thackeravlll970l) . Soft X -ray emission from 
the direction of HD 49798 was first detected with the 
Einstein Observatory in 1979, but it could be studied 
in detail only with a ROSAT satellite observation per- 
formed in 1992. The ROSAT data showed a very soft 
spectrum and led to the disc overy of a strong modula- 
tion with a period of 13.2 s ((Israel et al.l 119971 ). While 
the periodicity clearly indicated the presence of a com- 
pact object, the poorly constrained spectral fit resulted 
in a large uncertainty in the X-ray luminosity, making it 
impossible to distinguish between a neutron star and a 
white dwarf. 

Thanks to a long observation carried out with the 
XMM-Newton satellite in May 2008, we could establish 
that the subdwarfs companion is a white dwarf. Its 
X-ray luminosity, well constrained by the high quality 
XMM-Newton spectrum, is much smaller than that ex- 
pected from a neutron star accreting in the stellar wind 
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of HD 49798. The measurement of the X-ray pulse de- 
lays induced by the orbital motion, and the discovery of 
the X-ray eclipse, allowed us to derive the X-ray mass 
function and to constrain the system's inclination. This 
information, coupled to the accurately measured optical 
mass function, gives the masses of the two stars: M S( 2 = 
1.50 ± 0.05 M Q for the subdwarf HD 49798 and M WD 
= 1.28±0.05 M(t, for its white dw arf companion RX 
J0648.0-4418 ([Mereghetti et al.ll2009f) . The main param- 
eters of this binary system are summarized in Table [TJ 

RX J0648. 0-4418 is one of the most massive white 
dwarfs currently known and the one with the shortest 
spin period. This binary is also particularly interesting 
since systems of this kind, where a massive white dwarf 
accretes from a helium star, have been proposed as pos- 
sible progenitors of type la supernovae (jlben fc Tutukovl 
119941 : Wang et al. 2009). Here we report a comprehen- 
sive analysis of all the XMM-Newton observations and 
the results of a search for optical pulsations. 

2. X-RAY DATA ANALYSIS AND RESULTS 

XMM-Newton performed four short observations of 
RX J0648.0-4418 in 2002 and a longer one in 2008 (see 
Table [2]). To look for possible variability, the 2002 point- 
ings were taken at different orbital phasefl Based on 
the o ptical ephemeris of the system (jStickland fc Llovdl 
1994), the 2008 observation was scheduled to include the 
orbital phase 0.75, in order to look for the presence of 
an X-ray eclipse. This orbital phase was not covered by 
previous X-ray observations. 

Here we concentrate on data obtained with the EPIC 
instrument (0.1-12 ke V), consisting of two MPS and 
one pn CCD cameras (jTurner et al.ll2001t IStruder et al.1 
2001). During all the observations, the three cameras 
were operated in Full Frame mode (time resolution of 73 
ms and 2.6 s for pn and MOS, respectively) and with the 
medium optical blocking filter. The data were processed 
using SAS version 9.0. 

The source photons for the timing analysis were ex- 
tracted from a circular region of 30" radius and their ar- 
rival time corrected to the Solar System barycenter and 
for the orbital motion using the system parameters given 
in Table [TJ The source spectra were extracted from a 
slightly smaller circle (20" radius), in order to avoid pos- 
sible contamination at high energy from a source of sim- 
ilar brightness and harder spectrum located 70" away 
from RX J0648. 0-4418. The spectra were rebinned to 
have at least 30 counts per energy channel and to over- 
sample by a factor 3 the instrumental energy resolution. 
The background spectra were extracted from nearby re- 
gions, which were on the same chip as the source and 
where no sources were detected. Observation B was 
affected by high particle background, therefore we ex- 
cluded it from the spectral analysis. 

2.1. Variability 

The background subtracted light curves of RX 
J0648.0-4418 in the total (0.1-10 keV), soft (0.1-0.5 keV), 
and hard (0.5-10 keV) energy ranges are plotted as a 
function of the orbital phase in Fig. [TJ (pn) and Fig. [2] 
(sum of the two MOS). These figures clearly show that a 



significant flux is detected also during the X-ray eclipse, 
with net count rates of (2.8±0.3) x 10 -2 counts s^ 1 in the 
pn and (1.7±0.2) x 10 -2 counts s _1 in the sum of the two 
MOS. The presence of X-ray emission during the eclipse 
is confirmed by the image shown in Fig. [3[ which was 
accumulated selecting events in the time interval corre- 
sponding to the eclipse. 

Besides the variation due to the eclipse, the light curves 
shown in Figs. [TJ and [2] indicate the possible presence 
of some variability. In particular, the count rates in 
September 2002 (obs.D) were slightly larger than in the 
2008 observation, especially at energies above 0.5 keV. 
This variability is confirmed by the spectral analysis de- 
scribed below. 

2.2. Spectral analysis 

We first analyzed the X-ray emission detected dur- 
ing the eclipse, by extracting the pn and MOS spec- 
tra corresponding to the orbital phase interval 0.73- 
0.77. The three spectra were fitted together obtaining 
good results either with a power law with photon in- 
dex r=2.8±0.3 or with a thermal bremsstrahlung with 
temperature kT Br =0.55±gj keV (x 2 /dof=17.2/16 and 
24.7/16, respectively). In both cases we fixed the in- 
terstellar absorption to the value N#=2.7xl0 19 cm~ 2 
derived from the 2008 un-eclipsed pn data. The X-rays 
seen during the eclipse of the white dwarf are most likely 
due to emission from HD 49798, but even in the case of 
a different origin, we expect that this emission is present 
at all orbital phases. Therefore we included its contri- 
bution in all the subsequent fits to the non-eclipsed RX 
J0648.0-4418 emission. 

The 2008 out-of eclipse pn and MOS spectra were ex- 
tracted excluding the orbital phase interval 0.71-0.79 
and fitted jointly with different modelfl Single com- 
ponent models gave unacceptable results, while good fits 
could be obtained with two-component models consist- 
ing of a blackbody plus either a power law or a ther- 
mal bremsstrahlung. The best fit parameters are sum- 
marized in Table [3] The results differ slightly from 
those we reported earlier ()Mereghetti et al.l I2009D due 
to the inclusion in the present analysis of the eclipse 
component. The blackbody component, which accounts 
for most of the flux, is similar in both models, with 
temperature kTss^40 eV and bolometric luminosity 
Lbb ~8x10 31 erg s _1 (for d=650 pc). The hard compo- 
nent provides an additional luminosity of ~8xl0 30 erg 
s _1 (0.2-10 keV, corrected for the absorption). To check 
whether the uncertainties in the eclipse spectrum affect 
the overall results, we repeated the analysis with differ- 
ent bremsstrahlung temperatures of the fixed component 
in the range kTs r = 0.3-0.8 keV and found that the best 
fit values for the non-eclipse spectra change by less than 
10%. 

Also the spectra of the 2002 observations rule out single 
component models, but, owing to their lower statistics, 
the best fit parameters have larger uncertainties com- 
pared to those obtained from the 2008 data. To search 
for possible long term spectral variations we fitted the 
2002 spectra using a power-law plus blackbody model 



1 T he preliminary results of the 2002 observations ( |Tiengo et al.l 
2004) are superseded by those reported here. 



2 With the addition of the eclipse component with parameters 
fixed at the values of the bremsstrahlung best fit. 
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Fig. 1. — Background-subtracted orbital light curves of RX J0648. 0-4418 obtained with the pn camera in three energy ranges: 0.1-10 
keV (top), 0.1-0.5 keV (middle), and 0.5-10 keV (bottom). The labels indicate the five XMM-Newton observations listed in Tabled 
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Fig. 2. — Same as Fig.l, for the sum of the two MOS cameras. 
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with photon index, temperature, relative normalization, 
and absorption fixed at the 2008 best fit values. The 
eclipse component was always included, with all the pa- 
rameters fixed as described above. For simplicity we re- 
stricted the comparison to the pn data, which are those 
with the smallest statistical uncertainties. The results 
are shown in Figures [4] and [3 referring respectively to 
the observations of 2002 May (A+C) and 2002 Septem- 
ber (D). The residuals obtained in this way (middle panel 
of each figure) suggest that in September 2002 only the 
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power law component was brighter than in 2008, while in 
May 2002 the largest residuals are in the low-energy part 
of the spectrum, where the blackbody component domi- 
nates. Allowing the power law normalization to vary in 
the September 2002 fit, we obtained good results with a 
power law flux ^30% higher than in 2008 (see residuals 
in the lower panel of Fig|Sj). For the May 2002 spectrum, 
we kept all the power-law parameters fixed and allowed 
the blackbody component to change. This resulted in 
a best fit temperature kTss=31.6±0.2 eV and emitting 
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TABLE 1 

Main parameters of the HD 49798/RX J0648. 0-4418 binary system 



parameter reference^ 

Orbital Period P~^~ b 1.5476666±0.0000022 days (1) 

Eccentricity e (2) 

Inclination i 79-84° (1) 

Distance d 650±100 pc (3) 



HD 49798 

Mass M sd 1.50±0.05 M Q (1) 

Effective temperature T ef f 46,500 K (4) 

Magnitudes U=6.758, B=8.017, V=8.287 (5) 

Surface gravity log g 4.35 cgs (4) 

Radius R sd 1.45±0.25 R Q (3) 

Luminosity L 3xl0 37 ergs~ 1 (3) 



RX J0648. 0-4418 



Mass M WD 1.28±0.05 M Q 

Radius Rwd 3000 km 

Luminosity Lx 10 32 erg s _1 

Spin period P 13.18425±0. 00004 s 

Period derivative P -5xl0~ 13 s s" 1 < P < 9xl0~ 13 s i 



(1) 
(1) 



a References: (1) IMereghetti et"aTI lt2009T) ; (2) IStickland He Llovdl JT994T) ; (3) IKudritzki Ik Simonl (WfB) ; (4) IHamannl fZGffi l; (5) 
ILandolt He Uomotol (t2007lV 



TABLE 2 

Log of the X-ray observations of HD 49798/RX J0648. 0-4418. 



Obs. 


Date 


start-end 


Exposure 


Orbital 


Spin period 






(MJD) 


pn/MOS (ks) 


phase 


(s) 


A 


2002 May 03 


52397.46-52397.55 


4.5 / 7.2 


0.45-0.48 


13.18421(7)( a) 


B 


2002 May 04 


52397.98-52398.06 


1.4 / 5.6 


0.78-0.81 




C 


2002 May 04 


52398.56-52398.59 


0.6 / 2.5 


0.15-0.16 




D 


2002 Sep 17 


52534.58-52534.72 


6.9 / 11.9 


0.06-0.12 


13.1856(14) 


E 


2008 May 10-11 


54596.90-54597.38 


36.7 / 43.0 


0.56-0.87 


13.18425(4) 



( a ) Joint analysis of observations A, B, and C. 
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FlG. 3. — X-ray image (0.2-10 keV, pn camera) showing emission 
from HD 49798 during the orbital phase interval 0.73-0.77, corre- 
sponding to the white dwarf eclipse. The displayed region has a 
size of 1.5x1.5 arcmin 2 ; north to the top, east to the left. The 
circle, with a radius of 30", indicates HD 49798. The source to the 
SW of HD 49798 has a thermal X-ray spectrum and is probably 
associated to a late type star of magnitude V=13. 

radius Rbb=42.5±1.5 km (see residuals in lower panel of 
FigHJ. Table E] summarizes all the parameters obtained 
in this analysis of the 2002 observations. 

2.3. Timing analysis 

A best fit pulse period P = 13.184 25±0.00004 s was 
derive d for the 2008 observation by iMereghetti et al.1 
(2009). After correcting the times of arrival to the So- 
lar system barycenter and for the effects of the orbital 
motion of the system, we searched for pulsations in the 
2002 observations obtaining the values reported in Table 
O No long term variations in the spin period could be 
inferred by comparing the 2002 and 2008 values. A lin- 
ear fit gives a 90% c.l. interval -5xl0 -13 s s _1 < P < 
9xl0 -13 s s _1 for the period derivative. 

The folded light curves at soft and hard X-rays are 
shown in Fig. [51 No significant variations are seen by 
comparing the 2002 and 2008 lightcurves. On the other 
hand, all datasets show a clear difference between the two 
energy ranges. The soft band, where the blackbody com- 
ponent dominates, has a nearly sinusoidal profile with a 
pulsed fraction of 56%. Above 0.5 keV the light curve is 
instead characterized by two peaks, of unequal intensity 
and out of phase with respect to the soft pulse. 

2.4. Phase-resolved spectroscopy 

On the basis of the spectral results presented above, 
clearly showing the presence of two distinct components, 
we performed phase-resolved spectroscopy choosing dif- 
ferent phase intervals tailored to the light curves seen in 
the soft and hard energy ranges. We concentrate here 
on the 2008 data (excluding the eclipse). The 2002 data 
gave fully consistent, but less constrained, results, ow- 
ing to their smaller statistics. In all the phase-resolved 
spectral fits discussed below we included the fixed eclipse 
component, as was done for the phase-averaged analysis. 

We first searched for spectral variations as a function 
of the pulse phase in the soft component, by extracting 
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Fig. 4.— Spectrum of RX J0648.0-4418 obtained in May 2002 
(Obs. A+C) with the pn camera. Top panel: data and best fit 
power-law plus blackbody model (see parameters in Table[4jl. Mid- 
dle panel: residuals obtained with parameters fixed at the values 
of May 2008 (x 2 /dof = 61.9/32). Bottom panel: residuals of the 
best fit model (x 2 /dof = 41.8/32). 
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Fig. 5. — Spectrum of RX J0648. 0-4418 obtained in September 
2002 with the pn camera. Top panel: data and best fit power-law 
plus blackbody model (see parameters in Table [4{. Middle panel: 
residuals obtained with parameters fixed at the values of May 2008 
(x 2 /dof = 50.7/34). Bottom panel: residuals of the best fit model 
( X 2 /dof = 34.7/34). 

four spectra corresponding to two pairs of phase inter- 
vals selected on the basis of the 0.1-0.5 keV pulse profile: 
pulse minimum/maximum (phases 0.25-0.65/0.65-1.25), 
and pulse decline/rise (phases 0-0.4/0.4-1). All the spec- 
tra were fitted with the power-law plus blackbody model 
keeping N# linked to a common value. The blackbody 
temperatures obtained in the four spectra were all con- 
sistent with the value kTss~40 eV found in the phase- 
averaged spectrum. Acceptable \ 2 values were obtained 
by forcing a common value of kT b b and fitting together 
the minimum and maximum or the rise and decline spec- 
tra. Thus we conclude that all the variations in the low 
energy pulse profiles can be interpreted as a change in 
the intensity of the blackbody component that remains 
at a constant temperature. 

In order to investigate the phase dependence of the 
hard component, we extracted three spectra correspond- 
ing to the intervals marked in Fig. [5] first peak (PI, 
phases 0.15-0.35), second peak (P2, phases 0.65-0.85), 
and inter-peak (IP, phases 0.35-0.65 and 0.85-1.15). 
The three spectra were fitted with a blackbody plus 
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TABLE 3 

Results of phase-averaged spectroscopy of the 2008 observation (pn+MOS) 





Blackbody + Power law 


Blackbody + Brcmsstrahlung 


N H (cm-*) 


<2.4xl0 19 


<1.6xl0 19 


kT ss (eV) 




40.ol^;| 


(km) 


i7.9l?:? 


16.31^ 


r 


1.6±0.1 




F<5[ (erg cm- 2 s" 1 ) 


(1.7±0.1)xl0- 13 




kT Srem (keV) 




7.91^ 


F Brem ( er S Cm ~ 2 O 




(1.53±0.08)xl0- 13 


L<& (ergs-l) 


8.9xl0 31 


8.2xl0 31 


kT Bci (kcV) 


0.55< e > 


0.55( e > 


F<& (erg cm- 2 s- 1 ) 


4.3xl0- 14 < e > 


4.3xl0- 14 < e > 


X 2 / dof 


151.1 / 133 


158.7 / 133 



a Blackbody emission radius at infinity, for d=650 pc. 

b Observed flux of the power law component 0.2-10 keV. 

c Observed flux of the brcmsstrahlung component 0.2-10 keV. 

d Bolomctric luminosity of the blackbody component, for d=650 pc. 

e Fixed eclipse component 

All the errors arc at the 90% c.l. for a single interesting parameter 



TABLE 4 

Spectroscopy results of the 2002 observations (pn data) 



May 2002 (A+C) September 2002 (D) 



N H (cm- 2 ) 


2.7xl0 19 ( a ' 


2.7xl0 19 l a > 


kT BB (eV) 


31.6±0.2 


39.6( a > 


Rbb (km) 


42.5±1.5 


16.6< a > 


r 


1.57< a > 


1.57< a > 


(10- 13 erg cm- 2 s" 1 ) 


1.7(0) 


2.2±0.2 


Lbs (erg s -1 ) 


2.2 xlO 32 


8.4X10 31 


kT Bc! (keV) 


0.49< a > 


0.49< a > 


Feci (10- 14 erg cm" 2 s" 1 ) 


4.2< a ) 


4.2< a ) 


X 2 / dof 


41.8 / 31 


34.7 / 34 



a Fixed 

b Observed flux of the power law component 0.2-10 keV. 

All the errors are at the 90% c.l. for a single interesting parameter 
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power law model, keeping absorption and temperature 
fixed to the phase averaged values obtained with the pn 
(7V ff =2.7xl0 19 cm" 2 and kT BB =39.6 eV). The results 
are given in Table [5] The spectrum during the first peak 
is slightly softer than at other phases, as shown by the 
confidence level contours for the power law photon index 
and normalization plotted in FigjTJ 

3. SEARCH FOR OPTICAL PULSATIONS 

We observed HP 49798 /RXJ0648.0-44 18 with the 
Univ ersity of Cape Town CCD photometer ()0'Donoghuel 
1995) at the 1.9-m telescope of the South African Astro- 
nomical Observatory in Sutherland on December 25 and 
27, 2009. The observations were carried out in the U 
filter, with a time resolution of 4 and 6 seconds, respec- 
tively. The first observing run lasted 29.1 minutes, the 
second run lasted 60.2 minutes. 

A timing analysis of the optical data did not reveal 
any significant signal at, or close to, the frequency of 
the X-ray pulsations down to a level of 1.5 mmag (De- 
cember 25), and 0.6 mmag (December 27). These limits 
correspond respectively to a pulsed flux of 1.5x 10 -3 and 
0.6xl0~ 3 photons cm~ 2 s _1 A -1 at 3600 A, several or- 
ders of magnitude above the extrapolation to the optical 
U band of the X-ray blackbody spectrum. Even con- 
sidering the low energy extrapolation of the power law 
seen during Peak 1, i.e. the spectral component with the 
steepest slope (r=1.9±0.2), gives a 3600 A flux of a few 
10~ 6 photons cm' 2 s" 1 A' 1 , well below our upper limit. 
However, other mechanisms not directly related to the 
spectral components seen in X-rays might in principle 
contribute in the optical/UV band, and the above limit 
is the best currently available in the optical band for this 
pulsar. 



TABLE 5 

Results of phase-resolved spectroscopy of the 2008 
observation (pn data) 





Peak 1 


Peak 2 


Inter-peak 


N H (°) (cm" 2 ) 


2.7xl0 19 


2.7xl0 19 


2.7xl0 19 


kT ss M (eV) 


39.6 


39.6 


39.6 


R ss (km) 


15.1±0.8 


17.0±0.7 


16.6±0.4 


r 


1.9±0.2 


1.4±0.15 


1.39±0.15 


F p1 (10 -13 erg cm" 2 s" 1 ) 


1.55±0.25 


2.6±0.3 


1.2±0.2 


X 2 / dof 


38.4 / 28 


43.8 / 32 


83.7 / 62 



a Fixed 

b Observed flux of the power law component 0.2-10 keV. 

All the errors are at the 90% c.l. for a single interesting parameter 



4. THE X-RAY EMISSION FROM RX J0648.0-4418 

The results of X-ray observations with XMM-Newton 
spanning six years indicate that the properties of RX 
J0648. 0-4418 are rather stable. Only minor changes were 
seen in the spectrum, but it is not possible at this stage 
to attribute them to the different orbital phases of the 
observations or to a long term variability. The spectrum 
and luminosity seen with XMM-Newton are consistent 
with those measured with ROSAT in 1992 and Einstein 



in 1979, although the uncertainties in the latter data are 
rather large. No significant variations were seen in the 
spin period, pulsed fraction and shape of the pulse profile. 
Given that this is the only known X-ray emitting binary 
composed of a hot sub-dwarf and a white dwarf, it is 
worth examining different possibilities for the origin of 
the observed X-ray emission. 

The upper limit on the white dwarf secular spin-down 
(P <9xl0~ 13 s s -1 ), implies a rotational energy loss 
3 P ~ 1.5 xlO 36 erg s" 1 (for a moment of 



10 50 g cm 2 ). This is sufficiently high to 



l WD in 2 P 
inertia Iwd 

power the observed luminosity of ~10 32 (d/650 pc) 2 erg 
s' 1 , with an efficiency similar to that seen in rotation- 
powered neutron stars. Even if no convincing examples 
have been discovered, the existence of white dwarf ana- 
logues of rotation-power ed ne u tron stars is an intrigu- 
ing possibility (see, e.g., lUsovl (|1993[ )). Among isolated 
white dwarfs, even the one with the fastest ro t ationa l 
velocity (RE J0317-853, P=726 s, IVennes et ail (p003l) ) 
can develop a maximum potential drop along open field 
linefH of only AV = fl 2 B R 3 / c 2 -7xl0 9 (726 s/P) 2 
(B/660 MG) (R/5000 km) 3 V, despite its high magnetic 
field (B-170-660 MG). This should be compared with 
the much higher value ^10 15 -10 17 V typically reached 
in radio pulsars. White dwarfs in binary systems can 
rotate much more rapidly. After RX J0648. 0-4418, the 
next fastest rot ators are found in the dwarf nova WZ 
Sge (P=28.87 s, IPatterson et"all (11998ft) and in the i n- 
termediate polar AE Aqr (P=33 s, IPatterson! (| 1979ft ). 
Their magnetospheres can, in principle, accelerate par- 
ticles to high energies, but still few orders of magnitude 
below that reached in radio pulsars. A value AV ^10 14 
V could be attained in RX J0648.0-4418, but only in 
the case of a strong magnetic field, which, as we discuss 
below, seems rather unlikely. The possible detection of 
pulsed hard X-rays (10-30 keV) in AE Aqr has been 
interpreted as evidence for non-thermal p rocesses pow- 
ered by the white dwarf rotational energy (|Terada et al.l 
2008), implying an efficiency of ~0.1%. AE Aqr spins- 
down at a rate of 5.6xl0 -14 s s _1 , most likely as a re- 
sult of the magnetic propeller effect. Its K-type main 
sequence companion transfers mass trough Roche-lobe 
overflow, but the accretion stream is disrupted by the 
magnetic field of the white dwarf and most of the mass 
is ejected ()Wvnn et al.lll997D . AE Aqr exhibits a hard X- 
spectrum and strong rapid variability at all wavelengths, 
contrary to RX J0648. 0-4418 which shows a nearly con- 
stant flux, dominated by a thermal- like component. The 
significant differences, except for the short spin-period, 
between these systems disfavor the presence of a rotation- 
powered mechanisms in RX J0648. 0-4418. 

We believe that accretion provides a much more plau- 
sible and natural explanation for the X-ray emission ob- 
served from RX J0648.0-4418. The Roche lobe of HD 
49798 has a radius R L - a (0.38+0.2 Log (M sd /M WD )) 
— 3.2 i?Q, where a=5.5xl0 11 cm is the orbital sepa- 
ratio n. The sub-dwarf, with a radius R s( z=(l-45±0.25) 
R Q (|Kudritzki fc Simonl 11978ft . is significantly smaller 
than its Roche-lobe. Therefore, accretion can only oc- 
cur via stellar wind capture. Hot sub-dwarfs have rel- 

3 Assuming a dipole field 
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Fig. 6. — Folded light curves in the soft and hard energy ranges obtained in the observations of May 2002 (left panel), September 2002 
(middle panel) and May 2008 (right panel). 

flow onto the star's magnetic pole(s), and even prevent- 
ing the formation of an accretion disk in the systems 
with the strongest helds. Shocks are formed in the ac- 
cretion column above the white dwarf and the hot plasma 
in the post-shock region emits hard X-rays. This hard 
component is usuall y fit with mult i- tempe r ature plasma 
model s (see, e.g., [Cropper et al.l () 19990 ; lYuasa et al.l 
(HSU), but, as a first approximation, it can be described 
by a thermal bremsstrahlung with temperature kT~30 
(M/ Mq) Kg 1 keV, where M is the white dwarf mass 
and Rg its radius in units of 10 9 cm. Part of the hard X- 
ray emission may be intercepted and reprocessed by the 
white dwarf surface, giving rise to UV and soft X-ray 
thermal emission, with blackbody temperature of tens of 
eV. The relative importance of the soft and hard com- 
ponent depends on several factors, including geometrical 
effects, the white dwarf mass, magnetic field, accretion 
rate per unit surface, with the consequence that the ob- 
served values of the rat io L so f t / L^ard span a large range 
(|Evans & Heilierll2007[ ). 

The value L so f t / Lf lar d ~10 observed for RX J0648.0- 
4418 is rather high, but not uncommon. For example, 
similar or larger value s are found in EU UMa, D P Leo, 
and RX J1007-2016 (IRamsav fc Cropperl [200l . The 
blackbody temperature we derived for the soft compo- 
nent, as well as the dimensions of the emitting surface 
are well within the range of values measured for the po- 
lars. On the other hand, the best fit temperature of 
the bremsstrahlung component, kT tb ^8 keV is rather 
small if one considers that, based on the relation men- 
tioned above, massive white dwarfs should have the high- 
est shock temperatures. This discrepancy might have 
to do with the unique properties of this white dwarf, 
which accretes from a low density stellar wind rather 
than through Roche-lobe overflow. The intensity and 
configuration of the magnetic field is also relevant for 
the properties of the X-ray emission, and might explain 
the differences between this system and classical cata- 
clysmic variables. An upper limit on the dipolar field of 
RX J0648. 0-4418 can be set by the requirement that its 
magnetospheric radius, Ra/, be smaller than the coro- 
tation radius R CO R = (G M WD P 2 /4tt 2 ) 1 /3 = g x 10 8 



Fig. 7. — Confidence contour of the power law component pa- 
rameters in phase resolved spectroscopy. Npr, is the flux at 1 keV 
and r is the photon index. The curves correspond to the 68%, 90% 
and 99% confidence level. 

atively weak stellar winds, compared to main sequence 
or giant stars of the same spectral type, but HD 49798 
is one of the few sdO stars for which evidence of a sig- 
nificant mass loss has been reported. The P Cygni pro- 
files of its N V and C IV lines, indicate a mass loss in 
the range 5x 10 -10 -10 -8 M Q yr _1 and a wind termi- 
nal velocity v=1350 km s _1 , reached at only 1.7 R S( 2 
(|Hamann et al.lll981l ). A more precise estimate of M = 
3xl0 -9 Mq yr _1 has been obtained thanks to a recent 
analysis of th e optical /UV s pectra based on improved 
wind models (Hamann 2010]). The accretion rate onto 
the white dwarf can be estimated as Mwd^ M (R a /2a) 2 
= 8xKT 13 Mq yr™ 1 , where R a ~2GM TO /v 2 is the ac- 
cretion radius. The corresponding accretion luminosity 
is 3xl0 31 (Rwd/3000 km) -1 erg s _1 , in good agreement 
with the observed value. 

The presence of two components in the X-ray spec- 
trum of RX J0648. 0-4418 is similar to what is observed 
in many cataclysmic variables of the polar and interme- 
diate polar classes. In these systems, the magnetic field 
starts to dominate the motion of the accreting matter 
at large distances from the white dwarf, channeling the 
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cm. Since R« is denned by balancing the magnetic pres- 
sure and the ram pressure of the accreting matter, this 
translates in an upper limit on the white dwarf magnetic 
moment /i ^2xl0 29 G cm 3 . However, the above condi- 
tion does not exclude the possibility that the source be 
in the sub-sonic propeller regime, in which Rm^Rcqb , 
but the matter is too hot to accrete (|Davies et alJ fl979). 
Using the equation that defines the transition between 
the accretion regim e and the subsonic propeller regime 
(Bo zzo eEaD[2008), we obtain an upper limit /i<3 10 
G cm 3 , corresponding to a field of the order of 1 kG at 
the white dwarf's surface. 

5. HD 49798: THE FIRST X-RAY EMITTING HOT 
SUB-DWARF 

Early type stars emit X-rays with relatively soft ther- 
mal spectra (kT~0.5 keV) and luminosity proportional 
to the bolometric luminosity. The average relation Lx 
~10 -7 Lb i has been derived from the observations of 
main sequence, giant and supergiant stars of O spec- 
tral type. The X-ray emission is believed to originate 
from hot plasma heated by shocks a nd instabilities i n 
the strong stellar winds of these stars (|Pallavicinilll989l) . 
Despite their very high effective temperature, O type 
sub-dwarfs are much fainter than main sequence stars, 
and, up to now, X-ray observations provided only an up- 
per limit of ~10 31 erg s _1 for the sdO star BD -3°2179 
(jD'Antona et al.lll983l ). 

The X-rays we detected during the eclipse are well de- 
scribed with a bremsstrahlung of temperature 0.5 keV 
and luminosity 2xl0 30 erg s _1 . Given the bolometric lu- 
minosity of HD 49798 L6 o ;=10 39 L©, this corresponds to 
a ratio Lx /Lboi = 7x 10~ 8 . Thus, both from the spectral 
and luminosity point of view, the observed X-ray flux is 
consistent with emission from HD 49798. To our knowl- 
edge, this could be the first detection of X-ray emission 
from a sdO star. Further observations of this system, 
as well as searches for X-ray emission from other mem- 
bers of this class, can provide interesting information on 
the properties of X-ray production in stellar winds span- 
ning a large range of parameters, from hot sub-dwarfs to 
supergiants. 

6. DISCUSSION 

HD49798/RXJ0648.0-4418 is the only known X-ray 
binary containing a massive white dwarf accreting from 
a hot sub-dwarf. In principle, other systems of this kind, 
with a luminosity of ~10 32 erg s -1 , could be detected 
even at distances of several kpc by the current X-ray 
satellites. However, due to the very soft spectrum, the 
observed flux is very sensitive on the amount of inter- 
stellar absorption. For example, a column density of 
3xl0 21 cm" 2 would reduce by one order of magnitude 
the observed EPIC count rate of RX J0648.0-4418, mak- 
ing it difficult to identify other sources of this kind in 
the Galactic plane only through X-ray observations (RX 
J0648.0-4418 is at Galactic coordinates 1=254°, b=- 
19°). 

On the other hand, also from the point of view of the 
optical properties, this system seems quite uncommon. 
In fact, H D 49798 is ra ther atypical compared to other 
sdO stars (Heber 2009), which on average are less mas- 
sive (M^0.5 Mq), less luminous, and have a higher sur- 
face gravity. This is probably a result of its evolution in a 



close binary in which a common envelope phase occurred, 
as clearly indicated b y several arguments . The low H 
abundance (X w ~0.19. iKudritzki fe Simon! (|1978Q ). sug- 
gests that HD 49798 is the stripped core of an initially 
much more massive and larger star. Also the high abun- 
dance of N and low abundance of C confirm that its 
present surface layers once belonged to the outer part 
of the hydrogen-burning core of a massive star. Slightly 
different evolutionary scenarios have been proposed to 
explain this system. 

One possibility is that HD 49798 descend from a fairly 
massive progenitor that began to fill its Roche lobe be- 
fore helium ignition (Case B evolution) and is currently 
in a core He-b urning phase (jKudritzki fc Simonl 119781 : 
llben fc Tutukovl H993). In this case, the mass of the He 
star is related to that of its progenitor by M# e =0.043 
M 167 (jlben fc Tutukovl [19851 ) . implying that HD 49798 
was originally a star of ~8 Mq. Two other possibi l- 
ities have been considered by iBisscheroux et all (|1997|) . 
The first one is that HD 49798 consist of a degenerate CO 
core, surrounded by a helium envelope with a He-burning 
shell at its base. In this scenario, it would be the core of a 
progenitor of ^5 Mq which lost mass during a common 
envelope event when it was on the early asymptotic giant 
branch (AGB). Alternatively, HD 49798 could be burn- 
ing hydrogen in a shell, being the remnant of a star which 
shed its envelope while on the thermally pulsing AGB. 
However, this possibility is disfavored since it predicts a 
mass of only ~0.6 Mq for the sub-dwarf. 

At the observed mass transfer rate, with only a small 
(< 10~ 3 ) fraction of the total mass lost in the subd- 
warf's wind accreted, the white dwarf mass is increas- 
ing very slowly. At the end of the current He-burning 
phase, HD 49798 will expand and fill the Roche-lobe, 
giving rise to a much higher mass transfer rate. This will 
produce a higher X-ray luminosity and possibly bright 
outbu rsts like the one s een in the helium nova V445 Pup- 
pis (jWoudt et al.l 12009). which could be the descendent 
of a system similar to HD 49798/RX J0648.0-4418. The 
a mount of accreted mass , ~0.3-0.5 Mq (see figure 1 
of llben fe Tutukovl (|1994f )) . is sufficient to bring the, al- 
ready quite heavy, white dwarf above the Chandrasekhar 
limit. However, the fraction of mass that is effectively re- 
tained depends on the rate and composition of the accret- 
ing matter, on the mass, composition and temperature 
of the white dwarf, and on the poorly known relevance 
of other factors, such as, e.g., rotation, magnetic fields, 
wind outflows. Thus the fate of RX J0648.0-4418 is un- 
certain. 

Evolutionary computatio ns for RX J0648.0-4 418 have 
been recently performed (Wang fc Ha n] 120101 ) assum- 
ing the mass accumulation efficiency that takes into 
account the wind mass loss triggered by the He-shell 
flashes (jKato fc Hachisull2004| ). These indicate that RX 
J0648.0-4418 will reach a mass of 1.4 M© after only 
a few 10 4 years of Roche-lobe overflow, during which 
~(5-6)xl0 -6 Mq yr _1 of He-rich matter are steadily 
converted to C and O, while the unburned matter is 
ejected by th e syst em trough an optically thick wind. 
IWang fc Hani (|2010t ) assumed that RX J0648.0-4418 is 
a CO white dwarf that will explode as a type la su- 
pernova when the Chandrasekhar limit is reached. The 
fast rotation can increase the mass stability limit above 
the standard value for non-rotating stars, thus systems 
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like RX J0648. 0-4418 could be the progenitor of over- 
luminous type la supernovae. While an ONe composi- 
tion seems more likely for RX J0648.0-4418 in view of 
its large mass, also in this case there are some uncer- 
tainties and other factors that might play a crucial role. 
For example, it has been pointed out that massive white 
dwarfs c an have a CO composition as a result of rapid 
rotation (|Dominguez et al.lll993l I1996D . 

If RX J0648.0-4418 is an ONe white dwarf, an accre- 
tion induced collapse (AIC) might occur, instead than 
a type la supernova explosion, leading to the formation 
of a millisecond pulsar. There is strong evidence that 
MSP are old neutron stars which underwent magnetic 
field decay and spin-up due to accretion of mass and 
angular m omentum in low mass X-ray binary sy stems 
(LMXRB) (jBhattacharva fe van den Heuvejfl99lf ). This 
evolutionary scenario is supported by the observation of 
short period X-ray pulsations in the persistent and/or 
burst emission of several LMXRB, by the small eccentric- 
ity of MSP with white dwarf companions in the Galactic 
dislfl, and by the rece nt discovery of PSR J1023+0038 
([Archibald et al.ll2009l ). the long sought missing-link be- 
tween the LMXRB and MSP phase. On the other hand, 
the f ormation of MSP by AIC of accreting white dwarfs 
(e.g., iBailvn fc Grindlavl (|1990fl ) has been considered to 
reconcile the apparent discrepancy in the birthrate of 
LMXRB and MSP, and it cannot be excluded that a 
fraction of the MSP population is formed through AIC. 
The discovery of PSR J1903+0327, a Galactic plane 
MSP with a high eccentricity o rbit (e=0.44) around a 
solar mass main sequence star ([Champion et al.l 120081 : 
iFreire et al.ll20H . gives further evidence for the possi- 
bility of directly forming young MSP. The high spin rate 
and low magnetic field of RX J0648.0-4418 make it an 
ideal candidate for the direct formation of a MSP though 
AIC. 

7. CONCLUSIONS 

HD 49798/RX J0648.0-4418 is a unique system, but its 
X-ray properties are in several respects similar to those 
seen in white dwarfs of the polar and intermediate po- 
lar class. Its spectrum consists of a very soft, strongly 
pulsed, thermal like emission plus a harder component 
dominating above ~1 keV. The two components have 
different pulse profiles and show small uncorrelated vari- 
ations on long time scales. These similarities suggest 
that the X-ray emission properties depend mainly on the 
physical conditions near the white dwarf, with little influ- 
ence of the large scale accretion scenario (wind accretion 
in RX J0648. 0-4418 wrt Roche lobe overflow in cata- 
clysmic variables). 

With a dynamically measured mass of 1.28±0.05 Af Q , 
RX J0648. 0-4418 is one of the most massive white dwarf 
currently known and the one with the shortest spin pe- 
riod, only a factor ^5 larger than the break-up limit. It 
is also the only known white dwarf accreting from the 
wind of a hot sub-dwarf star. This system is the out- 
come of a common envelope evolution, most likely of an 

4 This is due to tidal circularization during the long LMXRB 
phase; eccentric MSP binaries can be formed if the companion is 
another neutron star that received a kick in the supernova explo- 
sion in which it was born, or for systems formed in dense envi- 
ronments where stellar dynamic interactions are important, like 
globular clusters. 



original pair of stars with mass of ~8-10 M Q . Its future 
evolution might lead to the formation of a neutron star 
through accretion induced collapse or to the explosion of 
a type la supernova. 

Both possibilities have interesting implications. Accre- 
tion induced collapse of a fast-spinning and low magnetic 
field white dwarf could be a promising scenario for the 
direct formation of non-recycled millisecond pulsars. In 
the second case, being the result of the evolution of rel- 
atively massive stars (~8-9 Mq), this would be a type 
la supernova formation channel with a short delay time. 
Future X-ray and optical observations can lead to more 
accurate determination of the system parameters, pro- 
viding crucial information to investigate in more detail 
the past and future evolution of this system. For exam- 
ple, further observations of the orbital phases including 
the eclipse can give a more precise measurement of its 
duration, a better characterization of the X-ray emission 
from HD 49798, as well as some information on the struc- 
ture of its wind by measuring possible variations of N# . 
High-resolution spectroscopy of the soft X-ray spectral 
component and deeper observations above a few keV are 
needed to better compare the properties of this peculiar 
X-ray binary with those of other accreting white dwarfs, 
and, in particular, to better constrain the magnetic field 
geometry and intensity. 
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